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ABSTRACT 

. ,— i | We estimate cosmological parameters using data obtained by the Very Small Array 

■ (VSA) in its extended configuration, in conjunction with a variety of other CMB 

' data and external priors. Within the flat ACDM model, we find that the inclusion of 

high resolution data from the VSA modifies the limits on the cosmological parameters 
as compared to those suggested by WMAP alone, while still remaining compatible 
with their estimates. We find that n h h 2 = 0.0234t°;°°^, n dm h 2 = O.lllj^gJe, h = 

0-73±g;g§, n s = 0.97lo; 3> 10 10 ^s = 23±^ and r = 0.14±g ^ for WMAP and VSA when 
no external prior is included. On extending the model to include a running spectral 
index of density fluctuations, we find that the inclusion of VSA data leads to a negative 
running at a level of more than 95% confidence (n run = —0.069 ± 0.032), something 
which is not significantly changed by the inclusion of a stringent prior on the Hubble 
constant. Inclusion of prior information from the 2dF galaxy redshift survey reduces 
the significance of the result by constraining the value of f2 m . We discuss the veracity 
of this result in the context of various systematic effects and also a broken spectral 
index model. We also constrain the fraction of neutrinos and find that f v < 0.087 at 
95% confidence which corresponds to m v < 0.32eV when all neutrino masses are the 
equal. Finally, we consider the global best fit within a general cosmological model with 
12 parameters and find consistency with other analyses available in the literature. The 
evidence for n mn < is only marginal within this model. 
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1 INTRODUCTION 

Recent measurements of the cosmic microwave background 
(CMB) anisotropies have allowed the determination of a 
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large number of cosmological parameters with unprece- 
dented accuracy. The Wilkinson Microwave Anisotropy 
Probe (WMAP) and pre-WMAP data sets can be fitted by 
a six-parameter ACDM model (see, for example, Bennett et 
al. 2003, Slosar et al. 2003). In order to break the degen- 
eracies inherent in the CMB power spectrum (Efstathiou & 
Bond 1999), various authors have augmented measurements 
of the CMB with observations of large-scale structure (LSS), 
for example, the 2dF galaxy power spectrum (Percival et al. 
2001, Percival et al. 2002), Lyman-a forest (Mandelbaum et 
al. 2003; Seljak et al. 2003), SDSS three-dimensional power 
spectrum (Tegmark et al. 2003), measurements of cosmic 
shear (Hoekstra et al. 2002) and the galaxy cluster luminos- 
ity function (Allen et al. 2003a), and/or information on the 
expansion rate of the Universe from measurements of the 
Hubble constant (Freedman et al. 2001) and high redshift 
supernovae (Perlmutter et al. 1999; Reiss et al. 1998). 

High-resolution (£ ^ 700) observations of CMB 
anisotropics provided by previously released data obtained 
by the Very Small Array (VSA, Grainge et al. 2003), the 
Arcminute Bolometer Array (ACBAR, Kuo et al. 2004) and 
the Cosmic Background Imager (CBI, Pearson et al. 2003), 
can also be important in reducing the impact of degeneracies 
and provide information on the parameters relating to the 
power spectrum of initial density fluctuations over a much 
wider range of scales. In particular, the WMAP team made 
use of these data in their analyses in order to improve the 
significance of their results (Spergel et al. 2003). 

In this paper, we study the cosmological implications 
of the new CMB power spectrum measured by the VSA 
which has a good signal-to-noise ratio out to a multipole of 
£ = 1500 (Dickinson et al. 2004). These observations cover 
33 fields, as opposed to 9 in Grainge et al. (2003), represent- 
ing an improvement of ~ 2 in signal-to-noise over the previ- 
ous data. By virtue of the accurately measured temperature 
of Jupiter by WMAP, the absolute calibration uncertainty 
for these data is reduced to 3% on the power spectrum; 
something which will be significant in our subsequent dis- 
cussion. The power spectrum is measured between £ = 300 
and £ = 1500 with a resolution in i?-space of A£ « 60. Previ- 
ous measurements of the power spectrum between £ = 130 
and £ = 900 using the VSA compact configuration can be 
found in Scott et al. (2003). 

We will first consider the standard six-parameter flat 
ACDM model, and then include extra parameters broadly 
in keeping with the approach taken in the papers published 
by the WMAP team (Spergel et al. 2003, Verde et al. 2003, 
Peiris et al. 2003). Our main focus will be on the initial spec- 
trum of fluctuations, quantified by the running of the spec- 
tral index, which appears to be particularly sensitive to high- 
resolution data such as ours. In the case where we do not 
impose external priors on the CMB data (WMAP+VSA), 
we find that there is significant evidence (> 2a) for negative 
running; something which is not implied by the WMAP data 
alone. The significance of this result is sensitive to the in- 
clusion of external priors, the relative calibration of WMAP 
and VSA, and possible source/cluster contamination of the 
measured power spectrum, illustrating issues which are of 
great relevance in the era of precision cosmology. The re- 
sult, if true, would be a significant challenge to models of 
slow-roll inflation, and so we also consider a broken spectral 
index model. As a final point, we consider a 12-parameter 



model fit to WMAP, WMAP+VSA and all available CMB 
data beyond £ > 1000, illustrating the effects of external 
priors on the estimated parameters. Our results within this 
model are compatible with previous determinations, both by 
the WMAP team and others. 



2 METHODOLOGY 
2.1 Cosmological model 

We will define the ACDM model as follows. First, we will 
assume that the Universe is flat and dominated by cold dark 
matter (CDM), baryons and a cosmological constant, A. The 
densities of these components relative to critical are denoted 
Odm, and £2a respectively and we define f2 m = Q dln + Q b 
to be the overall matter density (CDM and baryons) in the 
same units. The expansion rate is quantified in terms of 
the Hubble constant Ho — lOO/tkmsec -1 Mpc -1 and we al- 
low for instantaneous reionization at some epoch z rc (< 30) 
which can also be quantified in terms of an optical depth 
r. The so-called physical densities of the CDM and baryons 
are defined as uidm = Odm h 2 and uib = ilbh 2 . We will con- 
sider only adiabatic models and, guided by the predictions 
of slow-roll inflation, we parameterize the initial fluctuation 
spectrum of this model by 



P(k) = A s 



k 



(1) 



where k c — 0.05Mpc _1 is the arbitrarily chosen pivot point 
of the spectrum, ns is the spectral index and As is the scalar 
power spectrum normalization. 

We will modify this model by the inclusion of two other 
parameterizations of the power spectrum. We will, for the 
most part, consider a model with a running spectral index, 



''"•■)= Mr 



(2) 



so that the overall spectral index of fluctuations is a function 
of scale, ns(k), given by 



ns(k) 



rfflog P) 
d(log fc) 



ns + n r 



(3) 



where n run is known as the running of the spectral index. 
For slow roll inflation to be well defined, one requires that 
l^runl <S 1 1 — ns |/2 (Leach and Liddle, 2003). Under certain 
choices of priors we find that there is some evidence that 
this inequality is violated by the preferred fits to the data. 
We therefore consider an alternative model, which could 
be motivated by broken-scale invariance models of inflation 
(see, for example, Barriga et al. 2001), but is probably best 
thought of as a test of whether or not the data prefer a single 
power-law. The specific choice we will make is to consider 

n(k) — m for k < k c and n(k) — n,2 for k > k c , (4) 

with an appropriate normalization for k > k c so as to make 
the power spectrum continuous and the same value of k c as 
used in the standard ACDM model. 

In our discussion of systematic effects in section 4, we 
will consider the possibility of an extra component to the 
anisotropies with Ct = 2ttAx x 10~ 6 . Such a component is 
motivated by foreground effects due to point sources and the 
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Sunyaev-Zel'dovich (SZ) effect from galaxy clusters along 
the line of sight. The temperature anisotropics due to such 
a component will be (ATef = (l{l+l)Ci/(2n)) = A^i 2 and 
could be significant for t > 1000. The VSA has a sophisti- 
cated procedure to extract the effects of point sources using 
a dedicated, co-located, single-baseline interferometer (see 
Dickinson et al. (2004) for details), and the VSA fields have 
been chosen to avoid very luminous X-ray clusters. There 
could still, however, be some residual contamination. More- 
over, claims have been made of an excess signal between 
£ = 2000 and I = 4000 by the CBI team (Mason et al. 
2003), who attribute this to the SZ effect. If the signal is a 
large as is claimed, then it could be a contaminant even at 
lower £. By including such a component in the parameter 
fitting, it should be possible to constrain the contribution at 
I > 2000 as well as gaining some insights into the possible 
systematic effects of making such an error. 

The other parameters which we will consider in our 
analyses are: /„ = Q v /Qdm, the fraction of the dark mat- 
ter which is massive neutrinos; Ok = 1 — fitot (fitot = 
fidm + + f2„ + Oa), the curvature in units of the critical 
density; w = Pq/pq, the equation-of-state parameter for a 
dark energy component modelled as a slowly rolling scalar 
field; ut the spectral index of tensor fluctuations specified 
at the pivot point k c = 0.002 Mpc -1 ; R = At/As, the ratio 
of the amplitude of the scalar fluctuations, As, evaluated at 
k c — 0.05 Mpc -1 , and that of the tensor fluctuations evalu- 
ated at k c = 0.002 Mpc -1 . In addition to these parameters, 
for which we fit, we will also comment on various derived 
quantities: to, the age of the universe; erg, the amplitude of 
density fluctuations in the spheres of 8/i _1 Mpc. 

2.2 CMB data 

In this paper we will consider the cosmological implications 
of four different combinations of CMB data. 

• The first data set, denoted COBE+VSA contains the 
VSA data as described in Dickinson et al. (2004) com- 
bined with the COBE data (Smoot et al. 1992, Bennett 
et al. 1996). The purpose of this particular data set is to 
check the consistency of the VSA data with the concordant 
model, without imposing the strong constraining power of 
the WMAP data set (Bennett et al. 2003) . 

• The second data set, denoted WMAP contains only 
the WMAP temperature (TT) data (Hinshaw et al. 2003) 
and temperature-polarization cross-correlation (TE) data 
(Kogut et al. 2003). We use these data sets to provide a 
meaningful comparison with cosmological results obtained 
from other data sets, avoiding differences that might arise 
due to the priors and other methodological issues. 

• The third data set contains WMAP data and the new 
VSA data and is referred to as WMAP+VSA. In this data- 
set we supplement the accurate measurement of the first two 
peaks by the WMAP satellite with the VSA measurements 
of the power spectrum in the region between the third and 
fifth peaks. The importance of these data set is to illustrate 
the extra information that is available from the measure- 
ments of the power spectrum on small angular scales. 

• The last data set combines the previous two with all im- 
portant CMB experiments providing measurements in the 
region of the second peak of the spectrum and beyond, 



namely CBI, ACBAR, Boomerang, Maxima, DASI (Pearson 
et al. 2003; Kuo et al. 2004; Netterfield et al. 2002; Hanany 
et al. 2002 and Halverson et al. 2002, respectively). This last 
data set is hereafter referred to as AllCMB. 

Throughout our analysis we ignore small correlations 
between data sets that arise due to the fact that they have 
observed the same portions of the sky. This applies only to 
correlations between WMAP, that has used nearly all the 
sky, and terrestrial experiments, that have observed only 
small patches. In all cases, the decoupling of observed angu- 
lar scales and the fact that any given patch of sky observed 
by a terrestrial experiment makes up less than 1% of the 
WMAP sky coverage makes this approximation truly valid 
and far below systematic uncertainties. 

2.3 External priors 

In addition to the CMB data sets described above, we con- 
sider the effects of other cosmological data, not only to break 
the degeneracies, but also to see how the measured CMB 
power spectrum fits in the wider cosmological context. Each 
of these 'external priors' is discussed below. 

• The constraint on Hubble's constant obtained by im- 
posing the Hubble Space Telescope (HST) Key project value 
of H = 72 ± 8kmsec _1 Mpc -1 (Freedman et al. 2001) as a 
Gaussian distribution. The error-bar includes both statisti- 
cal and systematic uncertainty and prohibits the low density, 
low h universes allowed by the CMB data alone. 

• Constraints on large scale structure from the 2dF 
Galaxy Redshift Survey (Colless et al. 2001; Percival et al 
2001; Percival et al 2002), which provides measurements on 
scales 0.02 < fc/^tMpc -1 ) < 0.15. The 2dF data measure 
the power spectrum of the matter fluctuations in the lin- 
ear regime, which is linked to the spectrum of primordial 
fluctuations and the parameters of the standard model in 
a different manner to the CMB data and, thus, provide an 
important consistency check. 

• Constraints from Type la Supernovae (SNela) (Perl- 
mutter et al. 1999, Reiss et al. 1998), which help to break 
the CMB geometrical degeneracy and thus accurately deter- 
mine the ratio of matter to dark-energy components in our 
Universe. 

• Constraints from the gas fraction (f ga s) in dynamically 
relaxed clusters of galaxies (Allen et al. 2002) and from the 
observed local X-ray luminosity function (XLF) of galaxy 
clusters (Allen et al. 2003a). These data provide very accu- 
rate measurements of matter content of our universe, albeit 
with large systematic uncertainties. 

• Constraints from cosmic shear (CS) measurements 
(Hoekstra et al. 2002), which provide an independent re- 
striction in the £l m -ag plane from that implied by X-ray 
observations of clusters. 



2.4 Parameter estimation 

The parameter estimation has been performed using the 
COSMOMC computer package (Lewis & Bridle 2002) using 
the April 2003 version of the software (note that the de- 
fault parametrisation is different in the more recent versions 
of the COSMOMC package). The calculations were performed 
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Table 1. Priors used on each cosmological parameter when it is 
allowed to vary. The notation (a, b) for parameter x denotes a 
top-hat prior in the range a ^ x $C 6. 



Basic Parameter 


Prior 




(0.005,0.10) 




(0.01, 0.99) 
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on LAM clusters with a total of 42 CPUs at the IAC in La 
Laguna, Tenerife and the COSMOS supercomputer facility 
at the University of Cambridge. The COSMOMC software uses 
the Markov Chain Monte Carlo (MCMC) algorithm to ex- 
plore the hypercube of parameters on which we impose flat 
priors. These priors are listed in Table. 1. Additionally, the 
software automatically imposes the physical prior J7a > 0, 
which can significantly affect the marginalized probability 
distributions (see Slosar et al. 2003 for further discussion). 
For each considered model, we have run the software until 1 
in 25 of samples are accepted. Once this is achieved we ignore 
the first 200 accepted samples as a burn-in phase. In the flat 
models this leads to 65000 independent samples, and 200000 
in the non-flat case. These samples were then thinned by a 
factor 25 and used to plot marginalized probability distribu- 
tions with the GETDIST facility, which is part of the standard 
COSMOMC package. This program uses a smoothing kernel to 
infer a sufficiently smooth posterior probability curve from 
discrete MCMC samples. 

3 RESULTS 

3.1 Flat ACDM models 

3.1.1 Standard six-parameter model 

We begin our discussion in the context of the standard flat 
ACDM model with six free parameters (ui^, Ud m , h, ns, As, 
t), which was discussed in Spergel et al. (2003) for WMAP, 
with no external priors. We should note that it is, in fact, 
z TC which we allow to vary in our analysis, but we present r 
to be consistent with previous work. 

The marginalized distributions for the parameters are 
presented in Fig. 1 and the derived parameter estimates are 
tabulated in Table 2. The values for WMAP alone can be 
compared with those in Spergel et al. (2003). Noting that 
they present cu ln = f2 m /i 2 , instead of u>dm, there are only mi- 
nor discrepancies in the central values, although some of the 
limits appear to be somewhat larger. The preferred value of 
the redshift of reionization is z rc = 17lg. The inclusion of 
the high-resolution data from the VSA modifies the limits 




20 30 ffl 50 60 70 0,0 0.2 0.4 0.6 
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Figure 1. Marginalized distributions for the standard 6- 
parameter flat ACDM model with no external priors (that is, 
CMB alone) using COBE+VSA (solid-line), WMAP alone (dot- 
ted line) and WMAP+VSA (dashed line). 



on each of the parameters as one can see from Fig. 1 and 
these are most significant for ns , whose best fitting value re- 
duces from 1.00 to 0.97. The result for ns will be central to 
our subsequent discussion of the primordial power spectrum. 
The results from WMAP+VSA are very similar to those pre- 
sented in Spergel et al. (2003) for WMAP+ACBAR+CBI. 

We have also included in Fig. 1 the marginalized distri- 
butions and derived limits obtained from the COBE+VSA 
data set, all of which show compatibility with the re- 
sults of WMAP. One slightly unusual result is that for 
u)h which is much larger than the value suggested by 
WMAP, WMAP+VSA and standard Big Bang Nucleosyn- 
thesis, io h = 0.020 + 0.002, (Buries, Nollett and Turner 2001) 
and is a result of somewhat larger amplitude of the 3rd 
peak and the shifted first peak preferred by the VSA data 
(Rubino-Martin et al. 2003) in isolation (see Dickinson et al. 
2004 for a detailed discussion of the preferred peak structure 
of the current data). Comparing the derived distributions 
with those obtained by Slosar et al. 2003 (using the earlier 
VSA data presented in Grainge et al. 2003), we find that 
the results are fully consistent, but the additional VSA data 
has led to tighter parameter constraints. In particular, the 
upper limit on LO drn has been significantly reduced. 
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Table 2. Parameter estimates and 68% confidence limits for the 
standard six-parameter flat ACDM model. 



Parameter COBE+VSA 



WMAP 



WMAP+VSA 



h 

ns 

10 10 A S 



o.o328t8;88?? 

0.12518;8 3 2 ? 

77+ - 15 
u. 1 1 _ 17 

1 05+ ' 12 

251J 1 
Unconstrained 



0240+ - 0027 
u.u^4u_ 0016 

q 117 +0.018 
U.IK _ .018 

n 7 o+0.10 

1 00+ 09 

1 - uu _0.04 



2711 



n ,o+0.16 
U-l»_o. g 



n 0234+ - 0019 

U.UZJ4_q Q014 

q -,11+0.014 
u - ±1± -0.016 

n 7 o+0.09 

o.97l8;8 6 3 

o.i4t8£ 



Table 3. Limits on ns and n run in the flat ACDM model with a 
running spectral index for different CMB data sets and external 
priors. 



CMB 



External 



COBE+VSA 

WMAP 

WMAP+VSA 

COBE+VSA 

WMAP 

WMAP+VSA 

COBE+VSA 

WMAP 

WMAP+VSA 



None 
None 
None 

HST 
HST 
HST 

2dF 
2dF 
2dF 



0.93] 



0.94 



+0.07 
0.06 



0.96 



+0.07 
0.07 



-U.U81_ 04g 
-0 060+ 037 

U.UDU_g 036 

-0 069+ 032 
U.UD»_ Q 032 



0.92 



+0.11 
0.12 



-0.081 



r0.048 
-0.048 



n qK+0.06 



-0.060 



+0.037 
0.037 



0.93 



+0.06 
0.05 



1.00 



+0.12 
0.13 



0.95 



+0.05 
0.06 



0.93 



+0.05 
0.05 



-0 069+ 036 

U.UDM_ () q 36 



-o.o44t8:8S 

-o.038l8;8 2 3 f 
-0.049l8° 3 3 ! 



3.1.2 Running spectral index models 

In the previous section we saw that the inclusion of the VSA 
data to that of WMAP shifts the derived limits on the spec- 
tral index. Standard, slow-roll models of inflation predict 
that the spectral index will be a function of scale, albeit at 
a very low level, and it seems a sensible parameter to allow as 
the first beyond the standard model. The analysis of Spergel 
et al. (2003) and Peiris et al. (2003) provided evidence for 
a non-zero value of n Tun (— — 0.03ll8;8i?) when using CMB 
data from WMAP, ACBAR and CBI, along with large-scale 
structure data from the 2dF galaxy redshift survey and the 
Lyman-a forest. This result was discussed independently by 
Bridle et al. (2003), Barger et al. (2003), Leach & Liddle 
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Figure 2. Marginalized distributions for ns an d "run in the flat 
ACDM model with a running spectral index. Line-styles are as 
in Fig. 1. The external priors adopted are: none (top row), HST 
(middle row), 2dF (bottom row). 



(2003), Kinney et al. (2003), where it is was shown that it 
was highly dependent on the inclusion of the data from the 
Lyman-a forest, the veracity of which has been questioned 
(Seljak et al. 2003). 

We will start our discussion by considering the same 
model as in the previous section with no external priors, 
but with n rlm allowed to vary. The marginalized distribu- 
tions and derived limits on ns and n run are presented in the 
top row of Fig. 2 and the first three rows of Table 3 for 
COBE+VSA, WMAP and WMAP+VSA. The derived lim- 
its on Lo h , co dln and h are not changed appreciably and the 
other parameters, As and r (or z IC ) are strongly degenerate 
and 2 re will feature in our discussion below. 

The values of ns and not particularly well con- 

strained by COBE+VSA, but it is worth noting that even in 
this case there is a definite preference for a value of n rlm < 0. 
The results have been included for completeness and provide 
a useful cross-check. The results for WMAP are somewhat 
different to those presented in Spergel et al. (2003), some- 
thing to which we will return in the subsequent discussion. 
In particular we find that n run = — 0.060!q {j 3 ^' a 1-6<t pref- 
erence for n run < 0, as opposed to n rlm = —0.047 ± 0.04 
from Spergel et al. (2003). The significance of this result is 
improved to 2.2<r by the inclusion of the high resolution data 
from the VSA. These quantitative results are borne out on 
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Figure 3. As for Fig. 2, but for the parameters f2 m and re- 



examination of the likelihood curves. It is worth emphasizing 
that this result comes from CMB data alone. 

We have tested the sensitivity of this apparently strong 
result to the inclusion of external priors from the HST and 
2dF galaxy redshift survey, and the results are also presented 
in Fig. 2 and Table 3. We see that the effect of the HST prior 
is to relax marginally the constraint on n run , although there 
is a significant change in the derived limit on ns. We note 
that the results for WMAP alone are very similar with and 
without the HST prior. 

The inclusion of 2dF does significantly affect our results. 
Using just WMAP we find that there is only a marginal 
preference for n run < and the inclusion of VSA only yields 
a 1.4a result. We note that this is a shift in the derived 
value and the error bars do not change significantly; it is 
worth discussing the reason for this shift since it is due to 
the breaking of a degeneracy by the addition 2dF data. The 
main parameter combination which is constrained by the 
galaxy power spectrum is the shape parameter F — f2 m /t 
which arises from the size of the horizon at matter-radiation 
equality measured in redshift space. Hence, once combined 
with the CMB data the derived parameters f2 m and £Ia 
are constrained individually (Efstathiou et al. 2002). Fig. 3 
presents the marginalized distribution for these parameters 
for the three cases: no external prior, HST prior and 2dF 
prior. We see that for the first two cases, in which there is 
significant evidence for n run < 0, the preferred values of S7 m 
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Figure 4. Marginalized distributions for ni and 112 in the fiat 
ACDM model with a broken power-law index. The line styles arc 
as in Fig. 1. The left-hand column assumes the HST prior and 
the right-hand column assumed the 2dF prior. 

are much lower (extremely low in the no prior case) with 
the corresponding mean values of the distributions giving 
Qmh fa 0.17 — 0.18, whereas in the latter case f2 m « 0.3, 
h « 0.68 and fi m /i fa 0.21, closer to the value suggested 
by Percival et al. (2001) from the their analysis of the 2dF 
alone. 

We have also considered the effects of including other 
CMB information from the two other high resolution exper- 
iments ACBAR and CBI. We find that the inclusion of their 
results does not appear to be as significant as the VSA in 
preferring a value of n run < and that the result of consid- 
ering WMAP+ACBAR+CBI+VSA is very similar to just 
WMAP+VSA. We note that the ACBAR and CBI experi- 
ments quote large global calibration uncertainties (20% and 
10% in power), which we believe is at least as responsible 
for this result as their errors on the individual power spec- 
trum band powers. We note that the calibration uncertainty 
for the CBI is likely to improve for future data in a simi- 
lar way as for the VSA. A future re-calibration of CBI and 
more data are likely to address this issue (C. Contaldi, priv. 
comm.). 

3.1.3 Broken power law models 

In the previous section we have seen that there is some ev- 
idence for an initial power spectrum of density fluctuations 
which is not described by a single power law index. The 
running spectral index model is suggested by slow-roll infla- 
tion. However, the values which are preferred by the data, at 
least with some priors, are too large to come from standard 
slow-roll inflation and are incompatible with the idea of the 
spectral index being a power series in log(fc/fc c ). Here, we 
consider a model with two spectral indices n\ and ri2, with 
the cross-over point being k c — 0.05 Mpc -1 . 

The results obtained from this model are presented in 
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Fig. 4 for the HST and 2dF external priors respectively, 
using CMB data from WMAP and WMAP+VSA. We see 
that, in both cases and with both datasets, one obtains ni « 
1. The situation for ri2 is more complicated. For the HST 
case (left column) we see that the best fitting value is very 
low. In fact, it is lower than the lower limit we have included 
as a top hat prior. For WMAP values as large as n.2 = 1.4 
are not excluded, whereas the inclusion of the VSA has the 
effect of excluding models with n<i > 1. The inclusion of 
the 2dF prior (right column) has a strong effect, moving 
the distribution of n,2 to larger values, but still preferring 
n 2 < 1. In this case for WMAP we find that m = 0.99±0.04 
and ri2 = 0.97 ± 0.18 which suggests that something close 
to scale-invariant m = n,2 = 1 is preferred, whereas for 
WMAP+VSA m = 0.99 ±0.03 and n 2 = 0.88 ±0.15. While 
this is clearly compatible with scale invariant even at the la 
level, there is undoubtedly a preference for a broken spectral 
index when the VSA is included. 

It is clear that our results are compatible with those 
of the previous section on running spectral index models. 
Models with n Tun < have a lower value of the spectral 
index for k > k c than for k < k c and this is exactly what 
we find in this alternative parameterization. We should note 
that large variations in 712 only lead to much smaller changes 
in the actual power spectrum than one might expect naively 
from, for example, plotting the power spectrum for different 
values of ns within the standard ACDM model. 

3.1.4 Neutrino fraction 

As a final extension to our flat ACDM model, it is of inter- 
est to include the fraction /„ of dark matter in the form of 
neutrinos. Evidence for a neutrino oscillation, and hence for 
the existence of massive neutrinos, has been found by solar 
neutrino and atmospheric neutrino experiments (Fukuda et 
al. 1998, 2002; Allison et al. 1999; Ambrosio et al. 2000; Ah- 
mad et al. 2002). Further evidence for a non-zero value of the 
neutrino mass has recently been claimed from cosmological 
data (Allen et al. 2003b). 

In addition to obtaining constraints on /„ , the inclusion 
of this parameter will inevitably lead to some broadening of 
the marginalized distributions for the other parameters. Of 
particular interest is whether the constraints on the running 
spectral index derived above are robust to the inclusion of 
/„. We therefore include /„, with the top-hat prior given 
in Table 1, into the running spectral index model. In the 
analysis of this model, we include the 2dF external prior, 
since current CMB alone provide only a weak constraint on 

Fig. 5 shows the marginalized distributions obtained for 
f v , ns and n run for the three CMB data sets COBE±VSA, 
WMAP and WMAP±VSA. We find that the 95% upper 
limit provided by the COBE±VSA data set, /„ < 0.132, 
is only marginally larger than that obtained using WMAP 
data, f v < 0.090. The combination WMAP±VSA gives sim- 
ilar limits to WMAP, namely f v < 0.087, which corresponds 
to neutrino mass of ra„ < 0.32eV when the neutrino masses 
are degenerate. 

For the parameters ns and n Tun , we see that, as 
compared with those plotted in Fig. 2 (middle row), the 
marginalized distributions have indeed been shifted and 
broadened by the inclusion of /„ although the effects are 



not very strong. In particular, we note that our earlier find- 
ing of a preference for a non-zero value of n run has been 
weakened somewhat. A non-zero n run is still preferred, but 
at reduced significance. For the WMAP±VSA data set, we 
obtain n s = 0.94+° ° o l and n run = — 0.041±g;g|^ with 68% 
confidence limits. 

In the above analysis we used only 2dF as an external 
prior. It is of interest to investigate the effect of including dif- 
ferent combinations of the additional external priors listed 
in Table 1. The effect of these additional priors has been 
calculated by importance sampling our previous results. We 
also investigate the effect of including all recent CMB data 
into our analysis. In Fig. 6, we plot confidence limits on 
all the model parameters for each of our four CMB data 
sets, each of which, in turn, includes four different combi- 
nations of external priors: 2dF, 2dF±/ gas , 2dF±/ gas ±XLF, 
2df±HST and 2dF±CS. The points indicate the median of 
the corresponding marginalized distribution, and the error 
bars show the 68% central confidence limit. If the distribu- 
tion peaks at zero, the point is placed on the axis and the 
95% upper limit is shown. 

We see that the inclusion of the / gas and XLF external 
priors significantly reduces the error bars on all parameters. 
The most profound effect is obtained from the XLF prior for 
the parameters /„, <7s and z IC , as might be expected from 
Allen et al. (2003b). Indeed, it is only with the inclusion of 
the XLF prior that a non-zero value of /„ is preferred and 
only then at limited significance. For each of the CMB data 
set combinations, the best-fitting value in this case is /„ £s 
0.05, which corresponds to neutrino mass of m v « 0.18eV 
when the neutrino masses are degenerate, with a zero value 
excluded at around 96% confidence. For er 8 the inclusion of 
the XLF prior significantly reduces the best-fit value and 
the error bars for all CMB data set combinations. A similar, 
but less pronounced, effect is seen for z rc . 

3.2 General ACDM model 

Thus far we have considered only a limited range of flat 
ACDM models. In principle, one should properly include all 
the relevant unknowns into the analysis in order to obtain 
conservative confidence limits. In this section, we consider a 
more general ACDM model. In addition to including f v and 
n run , the standard six-parameter flat ACDM model is fur- 
ther extended by including Q. k , w, R = At/As and n-r- This 
gives 12 variable parameters in total, for which we adopt the 
top-hat priors listed in Table 1. 

For this model, we consider the three CMB data sets 
WMAP, WMAP±VSA and AllCMB. In addition, we now 
use both 2dF and SNela as our basic external priors, which 
are required in order to set constraints on our 12-dimensional 
cosmological parameter space. For each CMB data set, the 
marginalized distributions for each parameter are shown in 
Fig. 7. In addition, marginalized distributions are plotted for 
the derived parameters Qa, fim, to and as. The correspond- 
ing confidence limits on the parameter values are given in 
Table 4. 

In this more general model we see that the marginalized 
distributions of the parameters in our simpler models have 
broadened somewhat, but are still consistent with our earlier 
findings. Perhaps most interesting is the fact that some of 
the marginalized distributions change considerably as more 
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Figure 5. Marginalized distributions for /„, n s and n run in the extended fiat ACDM model, using the 2dF external prior and COBE+VSA 
(solid-line), WMAP alone (dotted line) and WMAP+VSA (dashed line). 
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Figure 6. Estimates for cosmological parameters in the flat ACDM running spectral index model, extended to include /„. Four CMB 
data sets are considered and, for each data set, four determinations are plotted, corresponding to different combinations of external 
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given by the HST key project, the Harrison-ZcFdovich value of the spectral index of fluctuations and a zero value for the running index. 



CMB data are used beyond WMAP. For Q,^h 2 we see a clear 
trend towards a lower preferred value (closer to the BBN es- 
timate) as one adds first VSA data and then all remaining 
CMB data sets. This effect is accompanied by a gradual 
upwards trend in the preferred Odm^ 2 value. The other pa- 
rameters exhibiting such trends are ns and n Iun . As more 
CMB data are included, the preferred value of ns moves 
slightly below unity, although this value is by no means ex- 
cluded. Perhaps more importantly, the upper limit on ns 
is significantly reduced as more CMB data are added. An 



analogous effect is observed for n run , for which the addition 
of VSA data significantly reduces the tail of the distribution 
for positive values of n run . 

The remaining marginalized distributions have very 
similar forms for each of the three CMB data sets, indi- 
cating that, for these parameters, the addition of further 
CMB beyond WMAP does not have a significant effect on 
their derived values and confidence limits. It is worth not- 
ing in passing, however, that all CMB data sets are fully 
consistent with a zero curvature model. Moreover, we find 
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Figure 7. Marginalized distributions for various cosmological parameters in the 12-parameter general non-flat ACDM model from 
WMAP (dotted line), WMAP+VSA (dashed line) and AllCMB (thick solid line) in combination with external priors from both 2dF 
and SNela. 



w = — 1 with an uncertainty of ± 24%, which is consistent 
with dark energy in the form of a cosmological constant. As 
regards inflation models, we find that the inclusion of VSA 
data makes a modest reduction in the upper limit on the 
tensor-to-scalar ratio, which is reduced still further (albeit 
marginally) by the inclusion of all CMB data; in this last 
case we obtain R < 0.68 at 95% confidence. The power-law 
index of tensor modes n T is fully consistent with zero. 

As we did for the fiat ACDM model, we may investi- 
gate the effect of including additional external priors in our 
analysis of the general model. In Fig. 8 we plot the confi- 
dence intervals on all the model parameters for each of our 
four CMB data sets, each of which, in turn, includes four 
different combinations of external priors: 2dF, 2dF+/g as , 
2dF+/ gas +XLF, 2dF+HST and 2dF+CS. Once again, we 
see that the inclusion of the / gas and XLF external priors 
has the greatest effect on the confidence limits, and that this 
is most pronounced for the XLF prior and the parameters 
/„, as and 2 rc . It is reassuring, however, that the derived 
limits on /„ for the general model are very similar to those 
obtained assuming the simpler flat model. We again find 
/„ ~ 0.05, with a zero-value excluded at about 92% con- 
fidence which is slightly lower than for the flat case. The 
effect of the XLF prior on as and z rc in the general model 
is also similar to that observed in the simpler flat case. 



4 DISCUSSION AND CONCLUSIONS 

We have used recent data from the Very Small Array, to- 
gether with other CMB datasets and external priors, to set 
constraints on cosmological parameters. We have considered 
both flat and non-flat ACDM models and the results are 
consistent. 

Within the flat ACDM model, we find that the inclusion 
of VSA data suggests that the initial fluctuation spectrum 
that is not described by a single power-law. As we have 
pointed out already, the value of n run preferred by the data 
is incompatible with the basic premises of slow-roll inflation. 
Moreover, the negative running, which reduces the amount 
of power on small scales and hence the amount of structure 
at early times, leads to predictions for the epoch of reioniza- 
tion at odds with the best fit to the CMB data. This comes 
almost directly from the temperature-polarization cross- 
correlation power spectrum observed by WMAP (Kogut et 
al. 2003). Given the implications of this result it is impor- 
tant to consider the possible systematic effects that might 
weaken it. 

The absolute calibration uncertainty of the VSA power 
spectrum is an important contributory factor to this result. 
The 3% uncertainty quoted in Dickinson et al. (2004) re- 
lies heavily on the the measurement of the temperature of 
Jupiter, Tjup, by WMAP and this requires an overall factor 
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of 0.92 modification in the power spectrum estimates from 
the previous VSA results (Grainge et al 2003; Scott et al 
2002) which were reliant on on earlier measurements of Tj up 
given by Mason et al. (1999). It was pointed out in Dickin- 
son et al. that, in fact, using an absolute calibration based 
on this measurement of Tj up gives the most consistent nor- 
malization of the power spectrum when compared to that of 
WMAP. 

We have investigated the effects of possible uncertain- 
ties in the calibration in two ways. First, we consider the 
possibility of using the Mason et al (1999) central value for 
Tjup while maintaining an overall uncertainty of 3%. As an 
alternative we just increase the overall uncertainty in the 
calibration to 10% while keeping the central value for Tj up 
from WMAP. The derived limits on ns and n run are pre- 
sented in Table 5 for these two possibilities using the HST 
and 2dF priors. We see in each case that the preference for 
n run < is weakened to below 2a compared to the calibra- 
tion based on WMAP's measurement of Tj up . It is clear that 
refinement of the absolute calibration of the VSA in the light 
of the WMAP measurements is something requiring further 
attention. 

Another possible systematic effect is the residual point 
source correction due to sources below out subtraction limit 



of 20m Jy. This was computed by normalizing the point 
source model of Toffolatti et al. (1998) to the observed VSA 
source counts which can then be extrapolated to lower flux 
densities. There are clearly some uncertainties in this proce- 
dure. It is possible that an imperfect subtraction, either an 
over-estimate or an under-estimate, could lead to inaccura- 
cies in the derived limits on the cosmological parameters, in 
particular on ns and n run . In order to investigate possible 
effects of such uncertainties we have performed our likeli- 
hood analysis with the inclusion of the parameter Ax which 
was discussed in section 2.1. We note that it is also possible 
that for Galactic foregrounds might contribute to this. How- 
ever, it was shown in Dickinson et al. (2004) that the level of 
foreground contamination of the VSA fields was negligible. 

We find that the derived limits on ns and n run in 
this case are less stringent than without including Ax for 
WMAP+VSA. The marginalized distributions for ns, n run 
and Ax are presented using CMB data from WMAP and 
WMAP+VSA for the external priors from HST and 2dF in 
Fig. 9 and the derived limits are presented in Table. 6. In 
fact the likelihood curves and derived limits for WMAP and 
WMAP+VSA are almost identical when Ax is included in 
the analysis and the WMAP limits are very similar to when 
Ax is constrained to be zero (see Table 3) . We see that there 
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Table 5. Limits on rig and n run in the flat ACDM model with a running spectral index for different absolute calibration schemes. The 
uncertainty refers to that in the power. See text for discussion. 
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Table 4. Parameter estimates and 68% confidence intervals for 
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Figure 9. Marginalized distributions for ns, ri run and Ax when 
Ax is allowed to vary. The line-styles are as in Fig. 1. but 
COBE+VSA is excluded The left-hand column is for the HST 
prior and the right-hand for the 2dF prior. 



are essentially no limits on Ax when just considering WMAP 
and that in the case of WMAP+VSA, A x is compatible 
with zero suggesting that to within at least w 100(/iK) 2 the 
source subtraction procedure has been successful. For the 
HST prior Ax < 214 (^K) 2 at 95% confidence and for the 
2dF prior A x < 155 (>K) 2 . 

Assuming that the subtraction is perfect and that the 
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Table 6. Limits on ns, "run and Ax in the flat ACDM model with a running spectral index when the parameter Ax is included. The 
final column is the 95% confidence upper limit on Ax- The units of Ax are in (pK) 2 . 



CMB External n s n Iun A x A x (2a) 



WMAP HST 0.95±gg| -0.059l^g^ Unconstrained 

WMAP+VSA HST 0.93t° -0.06lt°;°g -46 ± 132 < 214 



WMAP 2dF 0.96t5>;°|? -0.0361°;°^ Unconstrained 

WMAP+VSA 2dF 0.94+°;^ -0.043^°;°^ -86 ± 123 < 155 



SZ contribution to the power spectrum is oc £ 2 for £ < 4000 
then our results on Ax can be used to derive a limit on 
(ATe) 2 in a bandpower, B3000, covering 2000 «S £ 4000 as 
observed by CBI. Under the assumption that ATt oc £ 2 , we 
find that 

28 

B3000 = y Ax ■ (5) 

This leads to a limit of B3000 < 1997 (uK) 2 using the HST 
prior and B3000 < 1446 ("K) 2 for the 2dF prior. The value 
quoted by Mason et al. (2003) is B3000 = (508 ± 168) (uK) 2 
which is more stringent than our limit. However, this value 
and our limit are derived in very different ways. The mea- 
surement of Mason et al. (2003) is a direct limit from high 
resolution imaging of 3 deep fields. While it is direct, the 
measurement the global power spectrum could be signifi- 
cantly affected by sample variance; one could have observed 
a field which in there are more clusters than the global aver- 
age and hence obtain a biased estimation of the global power 
spectrum. Our limit is indirect, coming from the power spec- 
trum measured over 82 deg 2 at lower angular resolution and 
it requires that the power spectrum be oc £ 2 as well as the 
cosmological model to be correct. It is likely that this repre- 
sents a reliable upper bound since the power spectrum will 
grow less rapidly than £ 2 for £ > 1500 due to the fact that the 
clusters responsible for the SZ effect are not point sources. 
Moreover, it is not as sensitive to the Poisson distributed 
number of clusters in an individual field. A more realis- 
tic modelling of the SZ effect using accurate power spectra 
could yield a more stringent upper bound. The two methods 
provide useful complimentary information and it is possible 
that a much more stringent constraint on Ax will be pos- 
sible when the VSA observes with higher resolution in the 
near future. 

One important feature of the power spectrum observed 
by WMAP is the apparent absence of power at very low 
£. This could be due to some as yet unknown physics, or 
it could be a manifestation of the interaction between the 
subtle systematic effects caused by the side-lobes, the Galac- 
tic cut, and the power spectrum estimation algorithm used 
by the WMAP team (Efstathiou 2003). It is worth assess- 
ing to what extent our result is dependent on the measured 
anisotropies with £ < 10. By excluding the multipoles with 
£ < 10 from our analysis we find that ns = 1.01 ± 0.07 and 



n run = 0.007±0.049, strongly suggestive of an ns = 1, scale 
invariant initial power spectrum for WMAP with the 2dF 
prior, whereas n s = 0.97 ± 0.06 and n run = -0.015 ± 0.047 
for WMAP+VSA and the same prior. The weakening of the 
constraint on n run should not be a surprise since excluding 
multipoles with £ < 10 cuts out nearly a whole power of ten 
in k and n run is the coefficient of a power series in log(fc/fc c ). 
However, we see that the inclusion of the VSA tends to pre- 
fer a spectral index lower than just WMAP. It is clear from 
this that the reason for the preference for a negative run- 
ning of the power spectrum when multipoles with £ < 10 are 
included is the tension between the measurements at £ < 10 
by WMAP and for £ > 1000 by the VSA. 

We should comment briefly on one aspect of our analysis 
of the running spectral index models which is not ideal: the 
preferred values of z IC and r. Most recent analyses of CMB 
data include an upper bound one of these parameters. In 
Spergel et al. (2003) a flat prior of r < 0.3 was used when 
in some cases the data had a preference for a high value of 
t by virtue of the \ow-£ TE correlation power spectrum; our 
analysis is no different and we believe that this is responsible 
for the differences between our analysis and that of Spergel 
et al (2003). All the likelihood curves and derived limits 
have made the not unreasonable assumption that z rc < 30. 
However, in some cases, particularly those for which we have 
included no prior from 2dF, the preferred values of z IC are 
close to this limit, uncomfortably close in some cases and one 
might be concerned that our results are sensitive to this. 
We find that the models with n Tun significantly less than 
zero tend to have larger values of z TC which explains why 
our derived limit of n run = -0.060 ± 0.037 from WMAP is 
larger than the one quoted in Spergel et al (2003). It also 
suggests that, by excluding z TO > 30, we have weakened 
the constraint on n run rather than artificially modifying the 
preferred value away from zero. An epoch of reionization 
with z^ ~ 30 would seem unlikely in the context of early 
structures being the source of ionization, but it is clear that 
the data suggest it. 

For the general 12-parameter ACDM model, we find 
that our marginalized distributions for ns and n run are 
broadened, as one would expect. Nevertheless, even in this 
case, the addition of VSA data significantly reduces tails 
of the distributions for ns greater than unity and for posi- 
tive n run , as compared with using WMAP as the only CMB 
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data set. Indeed, these effects are reinforced by the use of the 
AllCMB data set. The inclusion of additional CMB data 
beyond WMAP also leads to a noticeable reduction in the 
preferred value of uj h and a corresponding increase in cj dm . 

To summarize, we find that there is evidence for n run < 
in a limited class of models, but within the general ACDM 
model with 12 parameters the evidence is much weaker. 
Standard models of inflation are generally incompatible with 
such large negative values of n run , but the data appears to 
point in that direction, although not totally conclusively. 
The inclusion of an external prior from 2dF appears to 
weaken the result by fixing f2 m w 0.3 in conjunction with 
the CMB data. The measurement of f2 m /i using the galaxy 
power spectrum is responsible for this shift. It is an interest- 
ing question as to how reliable this measurement is since a 
slight shift in the results, a preference for fl m h £s 0.17 rather 
than Q, m h « 0.21 would bring their preferred value into line 
with that suggested by the CMB alone and would uphold 
the possibility of n run < 0. Since none of the galaxy red- 
shift surveys have conclusively observed the turnover in the 
power spectrum on which this determination of Q, m h is based 
we assert that there is still room for some doubt. We have 
also investigated the possible systematic effects which could 
weaken our result. We believe that the assumptions behind 
the power spectrum measurements presented in Dickinson et 
al. (2004) are the best available using the observations which 
we have made and the other information from the literature 
we have used. For sure, we have shown that measurements 
of the CMB power spectrum beyond £ = 1000 can have an 
impact on the estimation of cosmological parameters and 
that future measurements in this region by the VSA, the 
PLANCK satellite and others will enable us in the future to 
make more definitive statements. 
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